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INTRODUCTION 

V-: / 

a — ^^A  computational  program  has  been  prepared  by  the 
Digital  Computing  Branch  to  provide  single  frequency  con- 
stant loss  contours  for  the  AMOS  near-surface  propagation 
model.  This  IBM  704  program,  designated  USL  Program  No. 
0131,  is  in  Fortran  language  and  Is  described  in  Appendixes 
A and  B. 
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THEORY 


Study  A of  reference  (a)  contains  a detailed  descrip- 
tion of  the  method  by  which  Project  AMOS  near-surface 
transmission  data  were  analyzed.  The  set  of  four  equa- 
tions which  were  developed  have  been  utilized  over  the 
past  decade  In  both  the  design  of  sonar  equipments  and  the 
development  of  optimum  tactical  procedures.  The  purpose 
of  Program  No.  0131  is  to  provide  an  Improved  method  for 
displaying  this  Information,  making  use  of  the  Laboratory'! 
automatic  X-Y  plotter  to  eliminate  the  need  for  laborious 
hand  graphing.  Contours  of  this  type  have  appeared  In 
Study  B of  reference  (a)  and  In  reference  (b).  It  Is  felt 
that  the  contents  of  reference  (a)  are  widely  enough  known 
in  the  sonar  field  to  preclude  repetition  here.  The 
remainder  of  this  memorandum,  therefore,  will  deal  with 
the  specific  details  of  the  program  design. 

Examination  of  the  equations  of  ref.  (a)  show  that 
most  of  the  parameters  are  functions  of  source,  receiver, 
and  layer  depth.  Therefore,  the  greatest  saving  In  com- 
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putlng  tine  can  be  nade  by  fixing  tbe  source  and  layer  depths 

I and  searching  along  a constant  receiver  depth  for  specific 

values  of  propagation  loss.  This  procedure  Is  simplified 
by  the  fact  that  the  loss,  for  these  reltr let  Ions , Is  Mono- 
tonia with  range.  The  only  complicating  factor  occurs  at 
two  zonal  boundaries,  at  which  a shift  'takes  place  from  one 
describing  equation  to  another.  Discontinuities  of  as  much 
as  several  db  have  been  observed  at  these  Interchange  points. 
To  avoid  possible  problems,  the  bounds  of  any  trial  range 
region  are  made  to  converge  until  they  sre  separated  by  a 
selectable  range  Increment , A R . Thereafter,  a linear  Inter- 
polation for  final  range  is  made,  regardless  of  the  total 
loss  variation  In  this  Interval.  The  program  accepts  values 
of  AR  between  0.1  and  1.0  kyds. 

The  limits  of  the  computation  are  set  by  establishing 
a "window”  In  range  and  depth,  which  corresponds  to  receiver 
locations  of  Interest  for  a fixed  source  depth.  Since  the 
AMOS  equations  assume  equal  loss  along  reverse  paths,  the 
nomenclature  "source  depth"  and  "receiver  depth"  are 
arbitrary,  and  the  output  curves  can  be  interpreted  as  show- 
ing the  effect  of  either  own  ship  or  target  motion  In  depth. 
The  range  Interval  of  the  window  is  bounded  by  "RmIN"  and 
"RMAX>"  where  RftfAX  BUlt  b*  no  greater  than  100  kyds,  and 
the  difference  no  less  than  10  kyds.  The  depth  Interval  Is 
bounded  by  "Z\(IN"  and  "ZjyiAX,"  which,  along  with  the  source, 
receiver,  and  layer  depths,  are  limited  to  a maximum  of 
1S00  feet.  The  depth  cuts  at  which  computations  are  made 
In  the  window  are  spaced  by  intervals  of  *AZ,"  which  can  be 
no  more  than  200  feet.  In  addition,  the  minimum  value  of 
4Z  Is  restricted  to  10  feet  or  0.01  of  the  total  window 
depth  Interval,  whichever  Is  larger.  The  set  of  constant 
loss  contours  are  bounded  by  "N^j^"  and  "Nmax".  wlth  inter- 
mediate values  spaced  at  Intervals  of  AN.  Twenty  such  con- 
tours may  be  computed  and  plotted  In  any  set;  however, 
only  the  first  thirteen  may  be  printed  out,  due  to  the  space 
limitation  of  the  printer.  Neither  N^]^  nor  may  exceed 

a value  of  200  db,  while  AN  Is  restricted  to  the  Interval 
between  1 and  20  db.  If  a contour  falls  outside  of  the 
window,  a "no  value"  Indication  Is  stored. 

The  remaining  parameters  which  affect  the  computation 
are  sea  state,  surface  temperature,  and  frequency.  Input 
temperatures  are  restricted  to  the  interval  between  30°  and 
100°F.  The  frequency  must  be  between  0.S  and  50  kc , even 
though  the  AMOS  data  represented  only  a 2 to  25  kc  variation. 


2 


It  1*  felt  that  these  extensions  are  somevhat  questionable, 
as  Is  the  maximum  depth  limit,  but  they  are  provided  as  a 
convenience  to  the  user,  assuming  that  he  Is  avare  of  the 
possible  errors  that  may  result. 

One  of  the  important  features  of  the  program  design  is 
an  automatic  checking  of  all  the  input  parameters  before 
the  computation  begins.  This  feature  prevents  the  loss  of 
computer  and  plotter  time  for  non-acceptable  cases,  and 
presents  the  user  with  a complete  listing  of  all  faulty 
Inputs.  An  example  of  such  a printout  Is  given  In  Figure 
(1)  . 

The  program  performs  all  computations  for  a single 
case  and  stores  the  resulting  range  points  In  a loss-depth 
matrix  prior  to  output.  Thereafter,  any  combination  of 
printing,  writing  tape,  and  punching  cards  can  be  selected 
by  the  settings  of  sense  switches  1,  3,  and  2,  respectively. 

COMPUTER  PROGRAM  DESCRIPTION 

A flow  chart  for  USL  Program  No.  0131  is  Appendix  A, 
and  the  IBM  704  Fortran  Program  Is  Appendix  B. 

The  basic  Input  data  deck  required  by  the  program 
consists  of  seven  cards  which  contain  the  Input  parameters. 
Consecutive  cases  can  be  run  by  adding  Indicator  cards  and 
either  one  or  seven  new  parameter  cards,  depending  upon 
the  extent  of  the  change  desired.  The  card  formats  are 
given  In  Table  1. 


TABLE  1 
CARD  FORMAT 


Card  No.  Cols.  Contents 


1 

1-4 

"0131" 

5-80 

(Unused ) 

2 

1-42 

Any  42  letter 

-number-symbol  combination 

43-80 

(Unused) 

3 

1-15 

rmin 

C<(Rmax  ~ 10 

kyds)] 

16-30 

rmax  <kyd8) 

C<100  kyds] 

31-45 

Afl  (kyds) 

C0.1  to  1.0 

kyds] 

46-80 

(Unused) 
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TABLE  l (Coni' d) 


CARD  FORMAT 


Card  No 

. Cols. 

Contents 

4 

1-13 

(ft)  [.<1500  ft] 

18-30 

(ft)  [<1300  ft] 

AZ  (ft)  [10  to  200  ft] 

31-43 

46-80 

(Unused ) 

3 

1-13 

Njujjj  (db)  [<200  db] 

16-30 

NmAX  <db)  [<200  db] 

31-48 

AN  (db)  [1  to  20  db] 

46-80 

(Unused ) 

6 

1-15 

S,  Sea  State 

16-30 

T,  Temperature  (°F.)  [30°  to  100°] 

31-80 

(Unused) 

7 

1-13 

f,  frequency  (ke)  [0.5  to  50  kc] 

16-30 

L,  layer  depth  (ft)  [^1500  ft] 

31-45 

Z0,  source  depth  (ft)  [«1500  ft] 

46-57 

Any  combination  of  letters  and  symbols 
used  to  Identify  card  output  packages. 

58-80 

(Unused) 
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(+)  means  that  cards  of  type  1-7  will 
be  read  In  for  the  next  case;  (-) 
means  that  only  a card  of  type  7 
follows . 

NOTE: 

All  numerical 

with 

the  exception  of  card  No.  1. 

The  program  use*  Tape  Unit  No.  3 for  primary  data 
Input,  and  Tape  Unit  No.  S for  output  to  the  plotter.  The 
functions  controlled  by  the  sense  switch  settings  are  listed 
at  the  beginning  of  the  Fortran  program  In  Appendix  B for 
convenient  reference.  At  least  one  of  the  first  three 
sense  switches  must  be  down,  or  no  output  at  all  will  result. 


If  a tap*  failure  ooours,  aa  option  la  provided,  through 
■ante  switch  0,  to  road  tha  Input  parameters  directly  from 
the  on-line  card  reader.  No  autonatlo  sealing  for  the 
plotter  is  Included  In  the  progran.  The  running  tine  Is 
dependent  on  the  else  of  AR  and  the  total  nunber  of  points 
computed,  whloh  Is  the  produet  of  the  loss  contour*  ( <20) 
and  depth  euta  (^100).  A typical  time  for  703  points  (13 
eontours  and  61  depths)  and  a AR  of  0.2  kyds  is  3 minutes 
for  tap*  only.  5 minutes  for  printout  and  tape,  and  60 
minute*  for  printout  and  punch.  The  us*  of  cards  as  either 
Input  or  output  node  should  be  discouraged,  since  the 
running  tine  per  case  becomes  obviously  prohibitive. 

Three  sample  computations  have  been  performed,  which 
correspond  to  the  parameters  of  Figs.  lA,  3A,  and  SA  of 
Study  B,  ref . (a).  The  printout  format  for  the  first  case 
Is  given  In  Fig.  2,  and  plots  of  the  three  cases  are  given 
In  Figs.  3,  4,  and  5,  respectively. 

SUMMARY 

An  IBM  704  Fortran  program,  USL  Program  No.  0131,  has 
been  written  to  provide  single  frequency  oonstant  loss 
contours  for  the  AMOS  near-surface  propagation  model.  This 
program  provides  an  Improved  capability  to  evaluate  the 
general  characteristics  of  an  extensive  sound  field  by 
automatically  plotting  the  oomputed  contours. 
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TESTING  DATA  REJECTION 
OUTPUT  TO  TAPE  NO,  5 


SOURCE  TOO  DEEP, 20  = 2225,0000000 


•AYER  TOO  DEEP.  L <=  2500.0000000 


FREQUENCY  ERROR,  F 


500.0000000 


TEMPERATURE  ERROR,  T * 212.0000000 


WINDOW  DEPTH  ERROR,  ZMIN 

ZMAX 


0. 

2000.0000000 


DELTA  Z TOO  SMALL,  DZ  * 


0. 


LOSS  DOMAIN  ERROR,  NMIN  =* 


NMAX 


DELTA  N IS  TOO  SMALL  AND 


0.5000000 


WINDOW  TOO  SMALL,  RMIN 


0.0010000  RMAX  - 


DELTA  R IS  TOO  LARGE  AND 


80.0000000 


CORRECT  THE  ABOVE  LISTED  INPUT  DATA  ERRORS 
THE  PROGRAM  CANNOT  8E  CONTINUED  FOR  THIS  DATA 

HIT  THE  START  KEY  IF  THE  DATA  FOR  A SUBSEQUENT  RUN  IS 

Fig.  1 - Printout  for  Automatic  Checking  of  Input  Data 


250.0000000 


0.2000000 


READY 
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SINGLE  FREQUENCY  AMOS-LOSS  CONTOURS  NO.  0131 
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APPENDIX  A 


FLOW  CHART  FOR  USL  PROGRAM  NO.  0131 
"CONSTANT  LOSS  CONTOURS  FOR  AMOS  PROPAGATION" 


12 


SUBROUTINE  "ANDY"  FOR  USL  PROGRAM  NO.  0131 


THIS  PI® 

<yfri 


IS  BBS1  Q.UALIW 
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APPENDIX  B 


FORTRAN  PROGRAM  NO.  0131 


SINGLE  FREQUENCY  AMOS-LOSS  CONTOURS  USING  SUBROUTINE  ANDY. 


SENSE  SWITCH  SETTINGS 


SENSE  SWITCH  1 (PRINTER) 
UP  - CONTOURS  NOT  PRINTED 


DOWN  - TABLE  OF  CONSTANT  AMOS-LOSS 
CONTOURS  IS  PRINTED 


SENSE  SWITCH  2 
UP  - NO  PUNCHING 


(CARDS) 


DOWN  - l.D.  CARD  AND  CONTOURS  ARE 
PUNCHED  FOR  DATA  PLOTTER 


SENSE  SWITCH  3 (OUTPUT  TAPE  UNIT  5) 


UP  - NO  OUTPUT 


DOWN  - CONTOURS  ARE  WRITTEN  ON  TAPE 
UNIT  5 FOR  THE  DATA  PLOTTER 


SENSE  SWITCH  A (NOT  USED) 


SENSE  SWITCH  5 


(NOT  USED) 


SENSE  SWITCH  6 (DATA  INPUT) 

UP  - INPUT  FROM  TAPE  UNIT  3 DOWN  - INPUT  FROM  CARD  READER  * 


* NOTE  THAT  THE  CARD  FORMAT  IS  IDENTICAL  IN  EITHER  INPUT  MOOE 


DIMENSION  RR(IOO) »EEN( 100) .R ( 100 , 20 ) ,B ( 20 ) 

EQUIVALENCE  (20.20) 

COMMON  R1.AA.G.H.AS.EL.2E.2E0.F.1N.20 
S ONE  DEC  1.0 
S TWO  DEC  0.000002 
S THR  DEC  3.0 
S ATE  DEC  8.000 
S TEN  DEC  0.001 
S 10D  DEC  131 
SVALN  OCT  000000777777 
SLIFT  DEC  80000 
SMONG  DEC  10000 
T 37  = 31.0 

I F ( SENSE  SWITCH  3)1000.100 
1000  REWIND  5 

END  FILE  5 
100  PRINT  1 

1 FORMAT  I 1H1 .30XA6HSINGLE  FREQUENCY  AMOS-LOSS  CONTOURS  NO.  0131) 
99  I F ( SENSE  SWITCH  6)98.97 
97  READ  INPUT  TAPE  3.96, ID 
96  FORMAT  I 14) 


1 1 »/)  «/)  l/)  i/l  i/)  i/) 


THIS  FAKE  IS  BEST  QUALITY  FJtACllCABL* 

J«OJ»  OQTY  VURUSHED  TO  DOC  __ 

S CLA  ID 

SUB  ! DO 
TZE*95 

294  PRINT  94 

94  FORMAT!  UH  WRONG  DATA  / 1H1  I 
HTRMOO 

95  READ  INPUT  TAPE  3 ,93 . HEAD , ER , CARD .NO .MEN ,CL A , TURE 
93  FORMAT (7A6) 

TRA*92 

98  READ  96. ID 

CLA  ID 
SUB  IDD 
TNZ*294 

READ  9 3. HE AD.ER.CARD.NO. MEN. CLA. TURE 
92  I F ( SENSE  SWITCH  3)201.202 

201  PRINT  206. HEAD. ER.CARD.NO, MEN. CLA, TURE 

206  FORMAT ( 1H0.32X7A6/1H  . 32X20HOUTPUT  TO  TAPE  NO.  5) 

GO  TO  204 

202  I F I SENSE  SWITCH  1)203.204 

203  PRINT  207, HEAD, ER.CARD.NO, MEN, CLA, TURE 

207  FORMA! ( 1H0.32X7A6) 

204  I F ( SENSE  SWITCH  2)205,89 

205  PUNCH  208, HEAD, ER.CARD.NO, MEN. CLA, TURE 

208  FORMAT (7A6) 

2 FORMAT ( 3F15. 7) 

89  I F I SENSE  SWITCH  6)86,87 

87  READ  INPUT  TAPE  3 . 2 , RM I N . RMAX , DELR .ZM I N ,ZMAX ,DELZ . ENM I N , ENMAX .DELE 
1N.SS.T 

85  READ  INPUT  TAPE  3 .83 , F ,EL .ZO.FL ,AG 
83  FORMAT  1 3F15.7  »2A6 ) 

82  SENSE  LIGHT  0 
10  = ABSF(ZO) 

IF(ZO  - 1500.0)500,500,501 

501  PRINT  502.  ZO 

502  FORMAT ( 23HOSOURCE  TOO  DEEP.ZO  = .F15.7/1H0I 
SENSE  LIGHT  1 

500  EL  = ABSFIEL) 

I F I EL  - 1.0)503.504,504 

503  EL  = 1.0 

504  I F ( EL  - 1500.0)505,505.506 

506  PRINT  507, EL 

507  FORMAT  ( 21H0LAYER  TOO  DEEP.  L =•  .F15.7/1HO) 

SENSE  LIGHT  1 

505  F = ABSF(F) 

IFIIF  - 50. 0)*<F  - 0.5001)508.508.509 

509  PRINT  510, F 

510  FORMAT ( 2 2H0FREQUENC Y ERROR,  F » , F15.7/1H0) 

SENSE  LIGHT  1 

508  SS  = ABSF(SS) 

T = ABSF(T) 

IF!  ( T-  30. 0)*(T-1  00. 0D511, 511. 512 

512  PRINT  513, T 

513  FORMAT! 24H0TtMPERATURE  ERROR,  T * .F15.7./1H0) 

SENSE  LIGHT  1 
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THIS  ?ASE  IS  BBS!  QUALITY  FRACIICABU 
FRO*  OOTY  JMRWSHKT  TO  DOC  - 


511  ZMlN  = ABSF(ZMIN) 

Z MAX  = ABSFIZMAX) 

IF(ZMIN  - 1500.0)514.514,515 

514  1F(ZMAX  - 1500.0)516.516.515 

515  PRINT  517  .ZM1N.ZMAX 

51/  FORMAT < 28H0W1NDOW  DEPTH  ERROR.  ZMIN  = .F15.7/20X8H  ZMAX  = .F15.7/1 

1H0) 

SENSE  LIGHT  1 

516  DELZ  = ABSF(DELZ) 

1FIDELZ  - MAX1F ( 10.0. ( IZMAX-ZMIN) /99.0 ) 1)518.519.519 

518  PRINT  520. DELZ 

520  FORMAT ( 25H0DLLTA  Z TOO  SMALL.  DZ  = .F15.7/1H0) 

SENSE  LIGHT  1 

519  I F ( DELZ  - 200.0)521,521.522 

522  DELZ  = 200.0 

521  ENMIN  = ABSF(ENMIN) 

ENMAX  = ABSF(ENMAX) 

IFIENMIN  - 200.0)523.523.524 

523  IFIENMAX  - 200.0)525,525.524 

524  PRINT  526. ENMIN. ENMAX 

526  FORMAT ( 27H0LOSS  DOMAIN  ERROR.  NMIN  « , F 1 5 . 7 . 5X6HNMAX  =.F15.7/1H0) 
SENSE  LIGHT  1 

525  DELEN=  ABSF(DELEN) 

IFIDELEN  -0.99998)527.527,528 

528  IFIDELEN  - 20. 0000 ) 529 . 5 30 . 530 

530  DtLEN  » 20.0 

529  GO  TO  531 

527  PRINT  532.DELEN 

532  F0RMATI27H0DELTA  N IS  TOO  SMALL  AND  *.F15.7/1H0> 

SENSE  LIGHT  1 

531  RMIN  = ABSF(RMIN) 

RMAX  = ABSF(RMAX) 

IFIRMIN  - 0.001)533,534.534 

533  RMIN  = 0.001 

534  1 F | RMAX  - 100.0)535.535.536 

536  RMAX  = 100.0 

5 35  I F ( ABSF (RMAX l-ABSF (RMIN)-IO.O I 537,538 .538 

537  PRINT  539, RMIN. RMAX 

539  FORMAT (25HOWI NDOW  TOO  SMALL.  RMIN  = , F 1 5 . 7 . 5X6HRMAX  =.F15,7/1H0) 
SENSE  LIGHT  1 

538  DELR  = ABSF(DELR) 

IF  I DEL R - 0.1)540.541.541 

540  DELR  = 0.1 

541  IFIDELR  - 1.0)542,542.543 

543  PRINT  544. DELR 

544  FORMA T ( 27H0DELT A R IS  TOO  LARGE  AND  =.F15.7/1H0) 

SENSE  LIGHT  1 

54?  I F ( SENSE  LIGHT  1)545.282 

545  PRINT  546 

546  FURMAI (43H0CORRECT  THE  ABOVE  LISTED  INPUT  DATA  ERRORS/46H0THE  PROG 
1 R AM  CANNOT  BE  CONTINUED  FOR  THIS  DAT  A / 1H0/60H0H I T THE  START  KEY  IF 
2 THE  DATA  FOR  A SUBSEQUENT  RUN  IS  READY/ 1H 1 ) 
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co  to  in  c/*  u>  o>  co  cn  in  cr>  co  co  co  co  Co 


THIS  TAUZ  IS  BEST  QUALITY  PRACTICABLI 
FHCX*  OCDTY  RtfttlSHBC  TO  GDC 


S HTR*81 

282  1 F ( SENSE  SWITCH  1)210.209 

210  PRINT  3.RM1N.RMAX.DELR.ZMIN.ZMAX.DELZ.ENMIN.ENMAX.DELEN.SS.T.F.EL. 

1ZO 

GO  TO  212 

209  I F ( SENSE  SWITCH  3)211.212 

211  PRINT  3.RMIN.RMAX.0ELR.ZMIN.ZMAX.DELZ.ENMI N.ENMAX .DELEN.SS.T.F.EL. 

1Z0 

212  I F ( SENSE  SWITCH  2)213.200 

213  PUNCH  208.FL.AG 
GO  TO  200 

86  READ  2 .RMIN .RMAX »DELR  *ZMI N « ZMAX .OELZ  » ENMI N * ENMAX .DELEN  * SS » T 
84  READ  83.F.EL.ZO.FL.AG 
S TRA*32 

81  I F ( SENSE  SWITCH  6)80.79 
80  READ  2. SIGN 
S T R A * 7 8 

79  READ  INPUT  TAPE  3.2.S1GN 
78  CLA  SIGN 
T PL  *99 

IFISENSE  SWITCH  6)84.85 
200  IZZ  = I ( ZMAX-ZM1 N ) /DELZ ) +1 .2 

INN  3 ( < ENMAX-ENMI N ) / DELEN )+l .2 
X = (6.0  - (2100. 0/{T  + 459.6)1) 

FT  * 1 • 23* ( 10. 0**X ) 

AA  = ( ( (0.651)*(F**2)*FT)  / ( ( F**2+FT**2 ) ) ) + I ( 0.02 69* ( F**2 ) > /FT ) 
ZEO*  SQRTF ( ZO/EL ) 

LOO  SS 
CLA  THR 
T LO*3  30 1 

AS  = 9.0*SQRTF (F/EL) 

GO  TO  3300 

301  AS  = 4. 5*SQRTF (F/EL ) 

3300  CONTINUE 
401  CONTINUE 
CLA  ZEO 
FSB  ONE 
TPL*12 

R 1 1 3 l } . O-ZEO) / 4 . 0 
TRA*13 

12  R 1 1 = SORTFI (ZEO**2 )-  1.01/5.0, 

13  Z = ZMIN 
101  ZE  = SQRTFIZ/EL) 

CLA  ZE 
FSB  ZEO 
SSP 

STO  TEMP 
FSB  ONE 
T ZE  * 1 A 
T PL  *1  A 

G = 0. 10*( 10. 0**(2.3*TEMP) )*( (F/25.0)**0. 3333333) 

S TRA* 1 5 

14  G = 20. 0*< (F/25.0)**0. 3333333) 

S 15  LOQ  T 37 
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THIS  PAGE  IS  BUST  QUALITY  mCTICABU 

CLA  Z E fUQ*  OOP!  JUOTISKED  TO  DDC  ' 

TLQ*411 
CLA  TEMP 
TLQ*411 
CLA  ZEO 
T LQ*4 1 1 

415  H = 0.4*( 10.0**TEMP+10.0**ZE+10.0**ZEO) 

TRA*4 1 2 

411  H = 4.0E+30 

412  CLA  F 
FSB  ATE 
TM I *2  20 1 

H = H*< (F/ATE)**0. 33333331 
2201  CLA  ZE 
FSB  ONE 
TPL* 1 6 

R1  = R11  + (1.0-ZEI/4.0 
S TRA*17 

16  R 1 = R 1 1 + SQRTF( (ZE**2)-1. 01/5.0 

17  RLBAR  * RM1N 

ENLBAR  = ANDYIZ. RLBAR) 

RUBAR  = RMAX 
ENUBAR  = ANDY1 Z, RUBAR > 
l = 0 
ENO  = ENM1N 
S 102  CLA  ENO 
S FSB  ENLBAR 

S TZE*1 8 

S STO  TEMP 

S SSP 

S FSB  TWO 

S T M I *1 8 

S CLA  TEMP 

S TM I * 1 03 

S CLA  ENO 

S FSB  ENUBAR 

S T ZE*19 

S STO  TEMP 

S SSP 

S FSB  TWO 

S TM  I * 1 9 

S CLA  TEMP 

S TM I * 1 05 

S CLA  1 

S TZE* 1 03 

S TM I *20 

RLBAR  = RUBAR 
ENLBAR  = ENUBAR 
RUBAR  =RR ( I ) 

ENUBAR  = EEN  < I ) 

I = 1-1 
S TRA* 1 02 

19  IZ  = ( ( Z - ZMINI/DELZ 1+1.2 

IN  = ((ENO  - ENMIN)/DELEN)+1.2 

R(IZ.IN)  = RUBAR 


20 


co  co  co  co  co  co  co 


S T R A* 1 04 

18  IZ  = ( (Z  - 
IN  = ((ENO 
R ( IZ . IN  ) = 
S TR A * 1 04 

103  IZ  = ( (Z  - 
IN  = ((ENO 
B R( IZ, IN)  = 

S 104  CLA  ENO 

S FAD  DELEN 

S STO  ENO 

S LDQ  ENMAX 

S T LQ* 1 06 

S T R A* 1 02 

105  IZ  = ( (Z  - 
IN  = ( ( ENO 
S CLA  RUBAR 

S FSB  RLBAR 

S FSB  DELR 

S STO  TEMP 

S SSP 

S FSB  TWO 

S TM I *2  1 


ZMINl/DELZ 1+1.2 

- ENMINl/OELENl+1.2 
RLBAR 

ZMINl/DELZ  1 + 1.2 

- ENMIN)/DELEN)+1.2 
VALN 


ZMINl/DELZ 1+1.2 
- ENMIN1/DELEN1+1.2 


S CLA  TEMP 

S TPL*22 

S TM I *2  0 

21  R ( I Z . I N 1 = RLBAR+ ( ( ENO-ENLBAR 1 / ( ENUBAR-ENLB AR ) 1 *DELR 
T RA* 1 04 

20  CLA  RUBAR 
FSB  RLBAR 
FSB  DELR 
LDQ  1 
HPR 

TRA*1 8 1 

22  ARE  = (RLBAR  + RUBAR1/2.0 
X * (ARE/DELR  + 0.5) 

CAY  = K 

ARE  = CAY*DELR 
EN  = ANDY  I Z .ARE ) 

CLA  EN 
FSB  ENO 
STO  TEMP 
SSP 

FSB  TWO 
TMI *24 
CLA  TEMP 
TMI *23 
TPL  *2  5 

23  RLBAR  = ARE 
ENLBAR=  EN 
GO  TO  105 

24  IZ  = <<Z  - ZMINl/DELZ 1+1.2 
IN  = ((ENO  - ENMINI/DELEN)+1.2 
R< IZ.IN)  = ARE 
RLBAR  * ARE 
ENLBAR  * EN 
GO  TO  104 
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(/'(/'CO  (/>  (/>(/»(/>  O (/»  CO  CO  (/>(/>  cn 


25  I = I + 1 

RR ( I ) = RUBAR 
EEN( 1 >*  ENUBAR 
RUBAR  = ARE 
ENUBAR*  EN 
GO  TO  105 
106  GLA  Z 

FAD  DELZ 
STO  Z 
LDQ  ZMAX 
TLQ*26 
T RA  * 1 0 1 

THIS  IS  LIKE  A SUBROUTINE 

26  I F < SENSE  SWITCH  11300.301 
301  I F I SENSE  SWITCH  3)302.303 
303  IFfSENSE  SWITCH  2)302.181 

181  DO  182  II  = 1,20 
STZ  BUI) 

DO  182  J=l,100 
STZ  R ( J , I I ) 

182  CONTINUE 
TRA*8 1 

300  Bill*  ENMIN 


JUIS  PAGE  IS  MSt  SJ^ITY  mCHCASL. 
yaofl  OOP*  WF0H3HB1  TO  DDC 


DO  30 A II  * 2. INN 
30A  BUI)  * B(II-l)  + DELEN 
PRINT  306 

PRINT  305,6(1) »B(2)»B(3),B(A),BI5),B(6),B(7), P'S), B(9).BIl0), B(ll) 
1.BI12I ,e(13) 

306  FORMAT ( 1HO.32X32HTA0UE  OF  CONTOUR  RANGES  IN  KYDS. ) 

305  FORMAT ( 1H0.6HZ  (FT) , F6 . 1 .2HDB . F6. 1 , 2HDB , F6. 1 ,2HDB • F6 . 1 » 2HDB , F6 . 1 ,2 
1H0B.F6.1,2H0B,F6.1,2HDB,F6.1,2HDB,F6.1,2HDB,F6.1,2HDB.F6.1,2HDB.F6 
2.1.2HDB.F6.1.2HDB) 

Z = ZMIN 
TRA*307  ' 

308  Z = Z + DELZ 

307  I Z = <<Z  - ZMIN) /DELZ  I +1.2 
IFUNN  - 13)3390.3390.3391 

3391  IPN  * 13 
GO  TO  3392 

3390  IPN  = INN 

3392  PRINT  309 ,Z . ( R ( I Z , II ) . II  * l.IPNl 

309  FORMAT! 1H0.F6. 1.13IF8. 2)  ) 

LDQ  Z 

CLA  ZMAX 
TLQ*308 
PRINT  1002 
1002  FORMA  T ( 1H 1 ) 


S T RA  *301 

302  EN  = ENMIN 
S TRA*3555 

28  EN  = EN  + DELEN 
3555  IFfSENSE  SWITCH  1)  31.35 
35  PRINT  135, EN 

135  FORMAT ( 1H0.3HN  =,F6.1»3H  DB) 
31  Z * ZMIN 
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THIS  PAGE  IS  BEST  QUALITY  PRACTICES  Li 
FBOa  OOPY  FURHSHKD  TO  GDC  — 

S TRA»3Z 

33  Z - Z + DELZ 

32  IZ  = (<Z  - ZMINJ/DELZ >+1.2 
IN  = ( ( EN-ENMIN ) /DELEN ) +1.2 
CAL  VALN 
COM 

ANA  R ( I Z * I N ) • 

TZE+344 

36  I F ( SENSE  SWITCH  2)37.38 

37  I P 1 = Z + 0.5 

IP2  = 100. 0*R ( IZ.INI+0.5 
I PI  = IP1  + MONO 
IP2  = I P2  + MONO 
PUNCH  901  . IP1 . IP2  .EN 

901  FORMAT (6XI6.6XI6.32XF15.7  ) 

38  I F ( SENSE  SWITCH  3)39.40 

39  I MX  = Z + 0.5 

I MX  = I MX  + MONG 
I MY  = 100  • 0*R ( IZ  . IN)  + 0.5 
I MY  = I MY  + MONO 

WRITE  OUTPUT  TAPE  5 . 1 004 . 1 MX . I MY 

1004  FORMAT (216) 

00  TO  40 

344  I F ( SENSE  SWITCH  1)40.34 

34  PRINT  11. Z 

40  CLA  ZMAX 
LOO  Z 
TLQ*3  3 

I F ( SENSE  SWITCH  3)1003.1006 
1003  WRITE  OUTPUT  TAPE  5 , 1 004 , L I FT » L I FT 

1006  I F ( SENSE  SWITCH  2)1005,1007 

1005  PUNCH  901  .LIFT, LIFT 

1007  CLA  ENMAX 
LOQ  EN 
TLQ*2  8 

I F ( SENSE  SWITCH  3)1881,181 
1881  END  FILE  5 
S TRA*1 8 1 

902  FORMA T ( 2F6 . 3 ) 

10  FORMAT ( 1H+.22XF15.7 .5XF15.7  ,/ ) 

11  FORMAT ( 1H  .22XF15.7.9X8HNO  VALUE./) 

3 FORMAT ( 9H0RIM1N)  =,F7.2.6H  ( K YD ) , 1 5X8HR ( MAX ) =.F7.2,6H  (KYDI.14X9H 
1DELTA  R =,F7.2,6H  I XYD) /9H0Z (MIN)  =,F8.2,5H  ( FT ) , 1 5X8HZ ( MAX ) =.F8. 
22.5H  (FT)  .14X9HDELTA  Z =,F8.2»5H  ( FT ) / 12HOLOSS (MIN)  =.F7.2.5H  (OB) 
3»13X11HL0SS(MAX)  =,F7.2»5H  ( OB ) » 12X1 3HLOSS  CHANGE  =.F7.2.5H  ( DS ) / 1 
42H0SEA  STATE  = , F 5 . 1 , 20X 13HT EMPER ATURE  =,F5.1.8H  (DEG  F ) / 1 2H0F REOUE 
5NCY  =»F5,1»5H  (KC) .15X13HLAYER  DEPTH  =.F7.1.5H  (FT), 10 X 14H SOURCE  0 
6EPTH  =»F7.1»5H  (FT) ) 

9 FORMAT ( 1H0.2XF15.7 ) 

ENDIO.l.  1.1,1) 
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FUNCTION  SUBPROGRAM  ANDY  0131 

FUNCTION  ANDY ( Z *R  1 
EQUIVALENCE  (ZO.ZO) 

COMMON  ARE1.AA.G.H.AS.EL.ZE.ZEO.F.IN.ZO 
SHALF  DEC  0.5000000 

S ONE  DEC  1.0000000 

S CON  DEC  0.43429448 

S ATE  DEC  8.0000000 

S SIX  DEC  60.0 
ZO  * ZO 

1 ARE  *=  R /SQRTF(EL) 

A3  TERM  = 20 . 0*CON*LOGF ( R ) + AA*R 
S CLA  ARE 

FSB  ARE  1 
TM I * 1 0 
FSB  HALF 
TZE*20 
TMI *20 
TRA*30 

10  R 1 = TERM  + 1 ARE/ARE1 ) *G  + 60.0 
S TRA*42 

20  R 1 = TERM  +12.0*1  ARE  - ARE1))*H  + 1 1 .0- 1 2. 0* I ARE  - ARE1)I)»G  +60.0 

S TRA*42 

30  R 1 = 1 1 TERM  - AA*R ) / 2 « 0 > + IAA+AS)*R  + H -AS* I ( ARE1+HALF ) * 1 EL** . 50  I 
1)  + 1 0 • 0*CON*LOGF 1 1 ARE1+HALF ) * 1 EL**. 50 ) ( + 60.0 

S 42  CLA  ONE 
S FSB  ZE 

S TMI *50 

S STO  El 

S CLA  ONE 

S FSB  ZEO 

S LRS  35 

S FMP  El 

S TMI *50 

ANDY  = R 1 
GO  TO  60 

50  El  = (25.0  -1 (ABSF1Z-EL) >**HALF)  - I 1 ABSF 1 ZO-EL ) ) **HALF I ♦ 5.0*R>* 
11  (F/25. 0**0.3333333) 

S CLA  El 

S TMI *912 

S 914  FAD  TERM 
S FAD  SIX 

S STO  El 

GO  TO  913 
S 912  STZ  El 
S CLA  El 

S T RA*9 1 4 

S 913  CLA  El 
S LDO  R 1 

S TLQ*51 
S STO  ANDY 

S T R A*60 

S 51  STO  ANDY 

60  RETURN 

END (0.1. 1.1.1) 
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